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ABSTRACT 
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Based on the MSS data taken in early Autumn by Landsat-1 and 2 of three 
years time lapse, Intensive studies on environmental change pattern in 
Central Japan Including Nagoya City, the 4th largest cfty of the country is 
made. , 

It is found that a patched cirrus reported by weather stations over the 
land which is hardly recognizable by eyes In the images of MSS affects' radiance 
values significantly In Band 4 while its effect is little In Bands 6 and 7 
(Near infrared spectra). The cross correlation coefficient analysis between^ 
two images indicates the highest value of 0.96 is obtained for a peninsula 
with. a smooth shoreline followed by a small island with the value of 0.95 
while the lowest value is 0.68 for a mountainous area with a river. The 
analysis of land use in Nagoya area shows there is little change in the 
metropol-i-tan area while fairly large change took place in the northern 
periphery of the city where large housing projects are going on. 


1. intruductron 

' Due to the fact that the radiance vl^ues obtained at the Landsat level 
is the energy of reflected solar radiation by earth's surface and scattered 
in the atmosphere, the effect of the atmosphere is doubled. Theoretical 
studies of atmospheric effects on solar radiation has been made by many 
investigators, among them are El terman (1 968) , Tanaka(197l) and Turner(l974) 
Rogers(1973) tried to evaluate atmospheric effect on the observed value of 
Landsat data, while Tsuchiya(l973, 1976) and Tschiya, Nakamura, Iwata and 
0chia?(1975) discussed the effect of thin cirrus on the radiance value of 
Landsat MSS data. 

Among a few good images taken over Japan during Lands 
is the one covering central Japan including Nagoya City 

Original photcgfapliy may be purchased froilis 
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From the view point of season this Image can be said that it was taken 20 
days before the Landsat-1 image which was taken on Oct. 5 1972 in spite of 
the fact that there is three years time lapse. Making use of these two images 
data an attempt is made to find out environmental change in Nagoya area. . 

The first subject discussed is the comparison of radiance values of four 
selected areas in reference with the meteorological conditions. The numerical 
analysis on the environmental change through computation of cross correlation 
coefficient between the two images described Before is the second subject- 
The third subject is the change of land use in Nagoya city and Its suburbs. 

i 

2. The image of central Japan taken by Landsat-1 and 2 

In Fig. ] (a - d) are shown the images of central Japan including 
Nagoya city, the fourth largest city of Japan with 2 million population. 

The images were taken by Landsat 1 and 2 on Oct. 5 1972 and Sep. il 1975 
respectively. Throughout the study the former is called the old while the 
latter is called the new for short. Comparing the new and old images a 
large differencejn the appearance is recognized over the ocean where' 
clouds are visible in the old image while not in the new. As far as a part 
of 'the~image over the land, no significant difference is recognized. 

3- Atmospheric condition 

In an attempt to see the actual atmospheric condition, weather map analy- 
sis is made and shown in Fig. 2, The area corresponding to the Landsat image 
is indicated by broken line. In the figure the weather, clouds and winds 
are plotted In standard WHO code. Numerals in the figure are visibility in 
the unit of km. 

On Oct. 5 1972 most of the stations reported patched thin cirrus and smog 
was observed at •Nagoya .. General cloud distribution was such that the clouds 
were thicker towards east, west and south. On Sep. 11 1975 no cloud was, 
observed at the weather stations within the image boundary. The distribution 
of visibility on both days is similar except eastern part where that of 
Oct. 5 1972 showed a poor visibility. The poor visibility of Nagoya is due 
to smog which is caused by dusts and smoke from the factories and heavy 
traffics. It is interesting to see that the visibility at Nagoya is same on 
both days, in order to see the vertical atmospheric condition, radio sounding 
at HanTamatsu (as to location see Fig. 2. (b)) are plotted and shown In Fig. 3. 
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Landsat MSS images taken on Oct. 5 1972 by Landsat-l and on 
Sept. 11 1975 by Landsat-2 respectively. The areas indicated 
by letters A, B, C and D are the areas over which digital 
analysis is made in section 4. 
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Fig, 3 Radio sounding at Haman 
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The vertical distribution of air temperature on both days Is similar 

except a very slight inversion in the lov/er layer close to the ground level 

observed on Oct. 5 1972, however there is a fairly large difference in 

the distribution of relative humidity which is shown in the right hand side 

of the figure. The -phase of vertical distribution is somehow opposite. 

There is an extremely dry layer at 800-mb level on September 11 1975 while 

on Oct. 5 1975 above 700-mb level it is very dry except at 500-mb. Qn 

Sep. 11 1975 a fairly thick wet layer exists between IfOO and 700-mb levels. 

/ ■ * - ' 

If aerosols in the upper layer is assumed same, the main factors contribut 
i ng to solar radiation reaching the ground is smog and water vapor. Since 
water vapor is a unique function of air timperature, it is considered that 
total amount of water vapor or precipitable water in the air column on both 
days is not much different in spite of the difference In vertical distribut- 
ion pattern of relative humidity. 

It. Radiance values of four selected areas- 

Making use of the CCT's supplied from NASA the mean • values and standard 
deviations of digital numbers corresponding to radiance over four selected 
areas, i . e. a farming area, Nagoya City, a mountain area with thick ever 
green trees and-sea surface of Ise Bay are computed and shown in Table 1. 

The locations of the areas are shown in Fig. 1 and 2. 

Table 1 Mean and standard deviation of digital number of four 

selected areas. A; Farming area, B: Nagoya City, C: Mountain 
area 'wi th"thlck Evergreen trees, D: Sea surface. 


- • 

. .00 GMT OCT 5 1972 

0048 GMT SEP 11 1975' 

B4 

B5 

B6 

B7 

B4 

B5 

B6 

B7 

A 

M 

S. D. 

39.3 

7.0 

32.0 

3.2 

43.5 

4.7 

21.5 

3.0 

22.8 

3.1 

22.9 

4.8 

47,1 

6.9 

21.7 

4.5 

B 

M 

S.D. 

37.8 

6.3 

31.7 

3.0 

29.1 

2.9 

12.1 

1.6 

27.7 

2.8 

30.5 

4.5 

31,1 

5.4 

10.9 
2.6 ‘ 

C 

M 

S.D. 

23.1 

8.3 

14.5 

3.2 

33.3 

6.7 

19.5 

4.5 

12.5 

2.0 

10.7 

2.9 

38.1 

8,4 

19.8 

5.0 

D 

M 

S.D. 

25.7 

1.4 

15.5 

1.3 

11.3 

1.4 

3*4 

0.8 

13.3 

0.8 

9.5 

0.8 

4.0 

0.9 

0.3 

0.0 
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Over the land areas, there is a certain similarity, i'. e. the digital 
numbers of Bands-4 and 5 of Oct. 5 1975 are systematically larger than those 
of Sep. 11 1975 even if the incoming solar radiation is a little smaller. On 
the other hand at Band 6 the value' of the latter is larger and at Band 7 the . 
value is nearly equal. Apparently the effect of thin cirrus and atmosphe'r ic 
condition are well reflected in the values of Bands-4 and 5 v/hile the contri- 
bution of them Is not so large in Bands-6 and 7 or it can be concluded that 
the effect is negligible. 

The global albedo defined as the ratio of outgoing and incoming solar 
radiation at the outside of the atmosphere is computed and shown in Table 
2. In the conversion of digital number into radiance, the values in Land- 
sat User's HandbooK( 1 972 ) is used. 


Table 2 Global albedo for the four selected areas {%) 



OCT. 5 1972 

SEPT. 11 1975 


B4 

B5. 

E6 

B7 

B4 

B5 

B6 

B7 

A' 

45-4 

36.0 

49;0 

49.1 

24.3 

23.9 

49.0 

45.5 

B 

43.6 

35.8 

32.9 

27.6 

29.3 

31,6 

32.2 

22.9 

C 

26.8 

16.4 

37.5 

44.5 

13.4 

11.1 

39.5 

41.6 

D 

l-2t0. 

17.5 

12. 6 

7.7 

14.2 

9.8 . 

. 4.1' 

0.0 


In average over the land the value of Band 4 of the old is larger than 
that of the new by I 6 % while the difference becomes -0.5 ?and 3-7 Z in 
Bands- 6 and 7 respectively. The larger value of the old in Bands- 6 and 7 
can be attributed to the effect of vegetation since the leaves of vegetation 
are more green with larger reflectance in near- infrared spectra. 

5 . Numerical analysis of environmental change during the three years 

In order to find out the degree of environmental change, a cross cor- 
relation analysis is made using the values of Band 7 data. Since the . 
radiance value has the least atmospheric effect as is indicated in the 
previous analysis. In the computation of cross correlation coefficient, 
the following procedures are taken. First, reference areas are selected 
from the old Image the location of which is indicated in Fig. 4. 
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were taken. 

From each reference area a template composed of 6^x64 picture elements 
(pixels) is made. Then cross correlation coefficients are computed between 
the template and the corresponding areas In the hew image. In order to 
avoid mis-registration, the computation Is made by Eq. (l). 

£ P [X(i+m, j +n) - X(m,n)][Y(i,j) - Y] 

C(m,n) “ ■■ - . — (1) 

\j w M ^ I M 1 

M E L [x(i+m, j+n)- -.X(m,n)]^ V ^ ^ - Y ] 

i=1 J-1 i-1 i=l 

where, '-X and Y denote radiance values of each pixel in the new and the old 
images respectively and M is fixed as 64. The size of the areas in the new' 
image Is fixed to be 4 times as large as a template with the template- 
approximately in the center. Then computation is made by changing 1 and m 
then the largest value of cross correlation coefficients is selected. The 
result of computation is indicated in Table 3. 
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The best fit is a peninsula 
with a smooth shore line (No 15) 
with the-value of 0.963 while 
the poorest value is 0.681j for 
a mountain area v/Ith a river. 

This poor value is considered 
due to the amount of water 
in the river since the amount 
of v/ater gives significant 
effect on the distribution 
of radiance especially in 
the near- infrared spectrum. 

As far as the numerical 
values are concerned, the areas 
with the values above 0.85 
may be defined as the area with 
little environmental change. It 
should be also added that the 
area named as a city (Mo 14) is 
picked up from eastern part of 
Nagoya city where land use change 


Table 3 The cross correlation coef- 
ficient between the old and the 
new i.mages for selected areas 


Number 

Geographical feature 

Cr. Cor. 
Coeff. 

1 

Harbor 

of Large city 

0.840 

2 

of Oil industry 

0.888 

3 

of Small city 

0.887 

4. 

Island 

Small 

0.928 

5 

Very small 

0.928 

6 

Large 

0,935 

- 7 

Small 

0,946 

8 

iMountain area with 

0.789 

9 

Ja small lake 

0.811 

10 

IMountain area with 

0.801 

11 

ja river 

0.684 

• 12 

Summit with trees 

0.809 

13 

Lake with a bridge 

0.935 

14 

Large city 

0,762 


Peninsula 


15 

with smooth shore 

0,963 


line 

16 

Same as above 

0:955 

17 

with Complicated 



shore ' line 

0.938 


took place during the three years. 


6. Land-use analysis in Nagoya City 

Further land use classification is made for Nagoya area. The area 

2 

analyzed- is approximately 15x15 km , Classification is made into 6 
categories, i, e. forest, agricultural, open, urban, suburban, water, new 
urban areas and others. The result of analysis is indicated in Fig. 5 
in which original 16 pixels are combined into one pixel therefore 
pixel In' the figure corresponds to the area of 316x316 km , In the figure 
a large area of blue color in the center is urban area and denotes the 
metropolitan area of Nagoya City. Comparing two images It can easily be 
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recognized Lhat very little change took place in metropolitan area while 
changes took place in the periphery of the city. 

To show more clearly the land use change, based on the two classified 
images, a land use change analysis is made and shown In Fig. 6 which shows 
large change in the areas to the north of metropolitan area. -Qrouhd surveys 
indicate that a number of large scale housing projects were going on at the 
observation time. it is also seen a fairly large variation of land use 
in the south-eastern area, however, here individual size Is not so large. 



Fig. 5 Land use classification of Nagoya City and its suburbs. 
Left hand: Oct. 5 1972, Right hand: Sep. 11 1975 
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Fig. 6 Land use change in Nagoya City and its suburbs during 
3 years between Oct. B 1972 and Sep. 11 1975. 
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7. Concluding Remark 

. The foregoing analyses lead to the following conclusion, in the digital 
analysis' of Landsat data, care must be taken for thin cirrus which often is 
hard to recognize visually. At the same time MSS data can be effectively 
used for the study of atmospheric radiation. Cross correlation analysis 
indicates the- .highest value obtained in central Japan is 0.963 for the area 
where little change occurred in land use. This suggests that a certa.in 
care must be taken in applying GCP matching method in precision processing 
of landsat raw data. Landsat MSS data are extremely useful for making land 
use map for a large area. 
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